Strong feedback and current noise in nanoelectromechanical systems 
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We demonstrate the feasibility of a strong feedback regime for a single-electron tunneling device 
weakly coupled to an underdamped single-mode oscillator. In this regime, mechanical oscillations 
are generated and the current is strongly modified whereas the current noise is parametrically big 
with respect to the Poisson value. This regime requires energy dependence of the tunnel amplitudes. 
For sufficiently fast tunnel rates the mechanical contribution to current noise can exceed the Poisson 
value even beyond the strong feedback regime. 
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Recent intensive research on nanoelectromechanical 
systems (NEMS) was motivated by a variety of physi- 
cal effects involved and the prospect of practical applica- 
tions 0. NEMS have been realized experimentally with 
molecules P|, semiconductor beams and suspended 
carbon nanotubes jj] . Phenomena observed include neg- 
ative differential resistance, phonon- assisted transport, 
and tuning the eigenmodes by the gate voltage. Most of 
these experiments were performed in the single electron 
tunnehng (SET) regime 0- 

In this regime, a NEMS is essentially a SET device 
coupled to a mechanical (harmonic) oscillator. The cou- 
pling is provided by a force F (see e.g. Ref. 0) acting 
on the oscillator, the value of the force depending on the 
charge state of the SET device. It determines the di- 
mensionless coupling parameter, A = F'^ /TiMuj^ , where 
M and lu are the mass and the frequency of the oscilla- 
tor. It was recognized already |3] that for strong coupling 
A ^ 1, mechanical degrees of freedom strongly influence 
transport through a SET device, leading, for instance, to 
polaron physics and Franck-Condon effect. However, the 
weak-coupling regime A <C 1 is characteristic for most of 
NEMS and will be considered below. 

Naively, the effect of the oscillator on transport cur- 
rent in this regime must be small and proportional to 
A. However, an underdamped oscillator can be swung 
up to big amplitudes even by a weak random force origi- 
nated from stochastic electron transfers through the de- 
vice 8] ; this amplitude may provide a strong feedback on 
the current. A less obvious effect is the extra dissipation 
due to electron tunneling which has been erroneously 
disregarded in We demonstrate in this Letter that 
such electron-induced dissipation may become negative, 
resulting in the generation of mechanical oscillation and 
in strong mechanical feedback. This takes place if the 
average charge accumulated in the SET device is a non- 
monotonous function of gate voltage. 

The strong feedback is the most manifest in the current 
noise. The natural measure of noise in nanostructures is 
the Poisson value 0, Sp = 2el. We demonstrate that 
in the strong feedback regime the noise is always para- 
metrically bigger than Sp due to long-time correlations of 



oscillator amplitude. If the generation is bistable, we pre- 
dict a telegraph noise that can be exponentially big. Even 
if the strong feedback is absent, the noise may still ex- 
ceed Sp. The experimental observation of the enhanced 
noise thus would provide a strong evidence for mechani- 
cal motion. 

SET systems are known to exhibit a (quasi) periodic 
structure of Coulomb diamonds in the plane of bias V 
and gate Vg voltages. Inside each diamond, the num- 
ber of extra electrons n is fixed to an integer [^. We 
concentrate on the region adjacent to the two neighbor- 
ing diamonds with n — and n — 1 where only these 
two charge states of the SET device participate in trans- 
port. In the classical limit, the statistical description of 
the system is provided by the joint distribution function 
Pn{x,v,t), with X and v being the position and veloc- 
ity of the oscillator. This distribution hmction obeys the 
following master equation (see e.g. Refs.8. .11,) . 
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Here, the total force F acting on the oscillator is the sum 
of the elastic force, friction force, and charge-dependent 
coupling force, respective to the order of terms in Eq. 
(0). Q ;» 1 is the quality factor. We count the position 
of the oscillator from its equilibrium position in the ri = 
state. In this case, F„ — nF. 

The "coUision integral" St [P] represents SET. There 
are four tunnel rates, P^^, where the subscripts L and 
R denote tunneling through the left or right junction, 
and the superscripts -I- and — correspond to the tun- 
neling to and from the island, respectively. In Eq. Ipj)l. 
P^ = P^ -|- P^. It is enough for our purposes to assume 
that each rate is a function of the corresponding energy 
cost AE^ associated with the addition/ removal of an 
electron to/from the island in the state n = 0/1 via left 
or right junction (Ai?^^ = —AEJ^p). Two independent 
energy differences are determined by electrostatics and 
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depend linearly on the voltages. Additionally, they are 
contributed by the shift of the oscillator, 

AE+ ^ ~W + Wl - Fx , AEj^ = -Wr + W + Fx , 

where we introduce a convenient parameter W represent- 
ing both eV and eV^, with Wl.Wr lying at the bound- 
aries of the diamonds and Wl <W< Wr in the trans- 
port region. The condition of applicability of the classical 
approach is that the energy differences are much bigger 
than energy quantum of the oscillator, W ^ huj. 

To simplify Eq. JQ), we implement the separation of 
the frequency scales: the inverse damping time k of the 
oscillator, the oscillator frequency w, and the total tun- 
neling rate = r+ + F^, assuming k ^ uj ^ Tt- The 
first condition implies that the mechanical energy hardly 
changes during an oscillation, while the second condition 
implies that the coordinate varies so slowly that F(a:) 
hardly changes between two successive tunneling events. 
In this case, we arrive at a Fokker-Planck equation for 
the distribution function of the slowest variable — me- 
chanical energy E, P{E). It reads 

'-^^CP,C^-^{E.iE)^DiE)^). (4) 

Here D{E) is the diffusion coefficient in energy space and 
the inverse damping time is given by k{E) = ii{E)+u}/Q, 
k being the SET contribution. It is instructive to express 
those parameters in terms of the average number of extra 
electrons in the island, n{x) = F+/Ft, 

r DiE)/E\_F^/J_( n(l-n)1\ 
\ k{E) j - M \Tt 1 dn/dW // ' 

Here, the angle brackets denote an average over the oscil- 
lation period, {A{x)) = J (de/n) cos'^ 9 A{x{E) sin 9), the 
oscillation amplitude being given by x{E) — ^J^eJm /oj. 

The SET contribution to the damping k has been er- 
roneously disregarded in Ref. |^ In fact, as Eq. © 
suggests, the diffusion and damping are closely related. 
In particular, in the absence of bias {Wl = Wr) the aver- 
age number of electrons is determined by the Boltzmann 
distribution and one proves that dn/dW = n{l — h)/kBT. 
In this case, the diffusion coefficient obeys the Einstein 
relation D{E) = kBTEn{E). This, in its own turn, guar- 
antees that Eq. Q is satisfied with the Boltzmann distri- 
bution P{E) oc ex:i>{-E/kBT). At eV > fcsT, the Ein- 
stein relation does not hold anymore. The effective tem- 
perature ED /k of the oscillator may become of the order 
of eV . Moreover, we will demonstrate that for energy- 
dependent tunneling rates, the damping k can become 
negative. This signals instability with respect to interac- 
tion with the oscillator. To stress the importance of the 
SET contribution we will disregard other contributions 
to the damping (Q — > oo), so that k = k. 



The stationary solution of Eq. Q apart from a nor- 
malization constant reads 

P(i;) cx cxp ^- ^ dE'E'K{E')/D{E')^ . (6) 

The current is modified by mechanical motion. At a 
given mechanical energy E, the current averaged over 
the oscillation period, Iw{E), is determined from the 
dependence of the current on the energy parameter 
I{W) = T+{W)T-^{W)/Tt{W) in the absence of oscil- 
lations: Iw{E) = J{d9/2T:)I{W + Fx{E)sm9). In the 
limit of small amplitudes, one has Iw{E) ~ I{W) + 
I"{W)XhujE/2. The actual current Iw is obtained by 
averaging Iw{E) over E with the distribution function 
P{E). Zero- frequency current noise in the Fokker-Planck 
framework is obtained as 

/•oo 

S=-4 dE5Iw{E)C-^5Iw{E)P{E), (7) 
Jo 

with 5I\\r{E) = I-w{E) — Iw ■ In our assumptions, the dis- 
tribution function is sharp at the energy scale of interest. 
Indeed, the typical mechanical energy needed to modify 
the rates is determined from the relation eV ~ Fx(E), 
yielding E ~ Miv'^^eV/ F)'^ . If A < 1, this always ex- 
ceeds the typical energy fluctuation eV. If the damping 
is positive at all E, P has a sharp maximum at E — and 
the current is very close to Iw{Q)- The average ampli- 
tude of the oscillations is too small to induce a noticeable 
mechanical feedback. 

The situation changes drastically if k{E) becomes neg- 
ative, indicating instability and growing amplitude of the 
oscillations. Since k{E) is determined by the tunnel rates 
only, the amplitude growth can only be stabilized by sig- 
nificant modification of the rates by the amplitude grow- 
ing: This is the strong mechanical feedback. Positions of 
probability maxima are determined by the roots of 

Ek.{E) = . (8) 

A non-trivial root Eq ^ indicates a generation of me- 
chanical oscillation with almost constant amplitude. This 
may strongly modify the current that is now given by 
IwiEo). Our analysis shows that the negative damping 
can only arise from the energy dependence of the tunnel 
amplitudes. This dependence is intrinsic for both semi- 
conductor quantum dots and molecules. 

To illustrate, we have chosen exponential energy de- 
pendence typical for wide tunnel barriers and one 
electron level in the SET system, 

F+^ - 2TlBe-''^-^<H^-M-AEl^)); 

= rO,^e'^-'-^^-"^«/;^(Ai?^_^) , (9) 

the factor 2 accounting for the spin degeneracy of the 
state n — I. The energy dependence sets a new energy 
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FIG. 1: Four stability regions in the gate-bias voltage plane. 
Bold solid lines indicate the edge of the Coulomb diamonds. 
Insets show the sketch of P{E) in each region. The horizontal 
lines indicate bias voltages used for current and noise scans 
in Figs. 2 and 3. 



scale Wc — i/ciL.R which is assumed to be much smaller 
than the charging energy. For concrete illustration, we 
choose UL = Om/Wc, ttR = 0.75/Wc, ksT = 0.2Wc, 
r" =^R- Fig. 1 presents the regions in gate-bias voltage 
plane that differ by number and stability of the roots of 
Eq. ((SJ. The region (i) corresponds to positive damping 
at all E and the absence of strong mechanical feedback. 
In the region (ii) the only stable solution corresponds to 
the generation of mechanical oscillations. There is bista- 
bility in the regions (iii) (stable roots at _E = and a 
finite amplitude) and (iv) (two stable roots at finite am- 
plitudes). Strong mechanical feedback is present in the 
regions (ii), (iii), and (iv). It is remarkable that the re- 
gion (iii) eventually extends to the Coulomb diamond 
where no current is possible at zero temperature with- 
out the oscillator: Generation of mechanical oscillation 
makes it possible. 

We illustrate the modification of the current by me- 
chanical motion in Fig. |21 The modification is noticeable 
provided the generation of oscillations takes place. It 
can be of the same order of magnitude as the unmodified 
current, and of either sign. The current can even exhibit 
jumps if there are at least two stable values Ei,E2 of 
the amplitude generated. The position of the jump cor- 
responds to the values of W at which the probabilities 
PiEi), P{E2) are equal. 

Let us now turn to the current noise. First we evaluate 
the noise in the absence of feedback (region (i) in Fig. 
1). We make use of Eq. (71) approximating D{E), k{E), 
5Iw{E) by their values at E 0. This yields 



/ d^iy D^jE) 
~ M^oj^ \dW^ J K^{E)E^ 



(10) 
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FIG. 2: Current modification in strong feedback regime for 
different bias voltages. The dashed (solid) lines give the cur- 
rent modified (unmodified) by mechanical motion. The two 
upper panels demonstrate that the modification is restricted 
to region (ii) where the generation takes place. The two lower 
panels illustrate the bistable regions (iii) and (iv). In the left 
lower panel, the current in the region (iv) gives a jump where 
the probabilities of two stable amplitudes are the same. To 
the right of the jump, the lower amplitude value is more prob- 
able. This value decreases and becomes zero at the border of 
region (iii). Therefore the modification of the current ceases 
there (see inset). In the right lower panel, the probabilities 
become the same in the region (iii). Therefore the modifica- 
tion ceases immediately after the jump. 



value Sp is of the order of 
Sp \ uo 



hujX 



(11) 



The small value of the second factor can be compensated 
by the large value of the first one. In this case, the electro- 
mechanical noise, concentrated at frequencies of the order 
of K, exceeds the Poisson value. 

In the region (ii), where the stable generation of the 
oscillation with Eq takes place, the current noise is due 
to small fluctuations of the oscillation amplitude. These 
fluctuations occur at a frequency scale of the order of 
k'{Eo)Eo- The noise is given by 



S' = 4 



I'HEo)DiEo) 



(12) 



_B=0 



The ratio of the electromechanical noise and the Poisson 



That is, it exceeds the Poisson value by a large factor. 
Our numerical results for regions (i) and (ii) (Fig. 3a, b) 
show that Eqs. 1)11(1 . (|12|l give a scale of the noise rather 
than a good estimation. The actual values of noise 
change by three-four orders of magnitude. The reason 
for that is that the parameters I' , I" , k, k' may become 
close to zero. 

In the regions (iii) and (iv) the oscillation amplitude 
randomly switches between two values Ei^2- Since they 
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2W/(eV) 



FIG. 3; Mechanical contribution to current noise for different 
bias voltages, (a) Noise in stable region (i) (Eq. IIUII becomes 
zero at /" = 0. (b) Noise in the region (ii) changes by orders 
of magnitude approaching zero at /' = and diverging at 
the boundaries of the region where Eo ^ (Eq. [Ell. (c,d) 
Telegraph noise is presented in the bistable regions only. The 
solid lines indicate the region where it is exponentially big. 
The cusp at maximum corresponds to equal probabilities of 
two stable amplitude values (As discussed, the current gives 
a jump at this point). 



correspond to two distinct values of the current /i_2, a 
telegraph noise is observed. The distribution function 
reaches maxima at E12 while the switching corresponds 
to passing the minimum of the distribution function, 
P{Em)- The switching times are therefore exponentially 
long. This may lead to exponentially big enhancement 
of noise. We compute the switching rates by Kramers' 
method 113, 



ri->2,2^i _ E,n>i'{E^) / 7(gi,2) P{E^) _ 

j{E)=EK'{E)/DiE), E^O; 
j{0)^2TrE^K^{E)/D{E)\E=o ■ 



(13) 



This yields the telegraph noise 



3(0) - 4(/i - 



(ri-2 + r2-i)3 



(14) 



that is exponentially big provided P{Ei)P{E2) > 
P'^{Em) (Fig. 3c,d). While this condition is not sat- 
isfied at the edges of bistability region, it certainly holds 
near the current jump, P{Ei) « P{E2), where the noise 
reaches maximum. The estimation reads 



To conclude, we analyzed the SET system weakly 
coupled to a mechanical oscillator and proved the ex- 
istence of significant modification of the current under 
condition of strong feedback where generation of me- 
chanical oscillation takes place. The latter is feasible 
for energy-dependent tunneling amplitudes. The current 
noise generated by mechanical motion in the strong feed- 
back regime significantly exceeds the Poisson value and 
may be exponentially big if the generation is bistable. 
Even if no generation takes place, this extra noise may 
exceed Sp for sufficiently fast tunneling rates. 
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